DISSIPATIVE NEUTRINO OSCILLATIONS 
IN RANDOMLY FLUCTUATING MATTER 



F. Benatti 

Dipartimento di Fisica Teorica, Universita di Trieste 
Strada Costiera 11, 34014 Trieste, Italy 
and 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste 



R. Floreanini 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste 
Dipartimento di Fisica Teorica, Universita di Trieste 
Strada Costiera 11, 34014 Trieste, Italy 



Abstract 



The generalized dynamics describing the propagation of neutrinos in 
randomly fluctuating media is analyzed: it takes into account matter- 
induced, decoherence phenomena that go beyond the standard MSW 
effect. A widely adopted density fluctuation pattern is found to be 
physically untenable: a more general model needs to be instead consid- 
ered, leading to flavor changing effective neutrino-matter interactions. 
They induce new, dissipative effects that modify the neutrino oscillation 
pattern in a way amenable to a direct experimental analysis. 
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1. INTRODUCTION 



When a neutrino propagates in a constant distribution of matter, coherent forward 
scattering phenomena can affect its time evolution. Despite the smaUness of the cross- 
section induced by the neutrino interaction with the medium, these matter effects can 
significantly modify the oscillation pattern, through the so-called MSW mechanism. [1, 2] 

However, forward scattering phenomena are just the simplest matter induced effects 
that can occur to a neutrino when the medium is allowed to fluctuate. In this case, the 
neutrino can be viewed as an open system, i. e. a subsystem immersed in an external envi- 
ronment (the medium); [3-7] its time evolution, obtained from the total ncutrino+matter 
dynamics by eliminating (i.e. integrating over) the matter degrees of freedom, is no longer 
unitary: it takes into account possible exchanges of entropy and energy between the neu- 
trino and the fluctuating medium. 

In many physical situations, one can safely ignore the details of the matter dynamics 
and use an effective description of the medium as a classical, random external field. Quite 
in general, any environment can be modeled in this way, provided the characteristic decay 
time of the associated correlations is sufficiently small with respect to the typical evolution 
time of the subsystem. In the case of relativistic neutrinos, this time scale can be roughly 
identified with the vacuum oscillation length: we shall therefore consider media that fiuc- 
tuate on time scales shorter than this. It has been recently pointed out that the interior 
of the Sim could indeed satisfy such a condition, [8] as likely as the earth mantel. Thus, a 
neutrino created in the sun or moving through the earth would effectively see a random 
fluctuating distribution of scattering centers and therefore be subjected to stochastic, in- 
coherent interaction with the medium. In this situation, correlations in the medium play 
a fundamental role: they are responsible for the generation of new matter effects, beyond 
the MSW ones, leading to irreversibility and loss of quantum coherence. 

The effects of fluctuating matter on neutrino propagation have been first discussed 
in [9-11], and recently reconsidered in [12-14]. However, all these analysis deal with a 
simple density fluctuation pattern, naturally suggested by the standard MSW treatment. 
Further, these fluctuations are assumed to be exactly (5-correlated; this is a highly idealized 
description of the environment, that e.g. for heat baths can be attained only in the limit 
of infinite temperature. 

Instead, in the following a more realistic exponentially damped form for the correla- 
tion functions in the medium will be adopted. Limiting for simplicity the discussion to the 
analysis of the oscillations of two species of neutrinos, we shall see that the effects induced 
by matter fluctuations can be fully described in terms of a limited number of phenomeno- 
logical parameters. They affect the oscillation pattern in a very distinctive way, that is 
amenable to a direct experimental study. 

On the other hand, when the simplified density fiuctuation hypothesis considered in 
[9-14] is adopted, a single constant is sufficient to parametrize the new matter effects. 
However, this approximation appears physically untenable, since by adopting it certain 
transition probabilities take unacceptable negative values; this serious inconsistency can 
be cured only by allowing more general matter fluctuations, pointing towards the presence 
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of flavor changing neutrino- matter interactions.^ 

As a final remark, it is interesting to point out that the dissipative effects induced by a 
randomly fiuctuating medium on neutrino oscillations involve in general the CP- violating 
phase that is present in the mixing matrix for Majorana neutrinos. Therefore, contrary to 
the vacuum case, matter oscillation experiments can provide, at least in principle, a way 
to distinguish between Dirac and Majorana neutrinos. 



2. MASTER EQUATION 

In discussing the mixing of two neutrino species, we shall adopt the familiar effective 
description in terms of a two-dimensional Hilbert space; [18-22] the fiavor states, that we 
shall conventionally call and It'p), will be chosen as basis states. With respect to 
this basis, the physical neutrino states are then represented by density matrices i?, i.e. 
by hermitian 2x2 matrices, with non-negative eigenvalues and unit trace. Their time 
evolution equation can be cast in a standard Liouville - von Neumann form: [23, 24] 

^ = -i[Ho.R{t)\+U[R{t)] . (2.1) 

The first piece on the r.h.s. describes the propagation of the neutrinos in vacuum; in the 
chosen basis, the effective hamiltonian Hq takes the standard form: 

Ho = ajfi-a, (2.2) 

where uj = Am^/4£', Am? being the square mass difference of the two mass eigenstates 
and E the average neutrino energy, while the unit vector n = (sin 2^, 0, — cos 2^) contains 
the dependence on the mixing angle, a = (o"i, cr2, (T3) being the vector of Pauli matrices. 
The additional contribution Lt[R] takes into account the presence of matter. As explained 
above, we shall consider the case of a rapidly fiuctuating medium, which can be described 
by classical stochastic fields. Its action on the travelling neutrinos can then be expressed 
via the commutator with a time-dependent hermitian matrix V{t), 

Lt [R{t)] = -i [V{t), R{t)] , V{t) = V{t) ■ a , (2.3) 

whose components Vi{t), V2(t), V^it) form a real, stationary Gaussian stochastic field V{t)] 
they are assumed to have in general a nonzero constant mean and translationally invariant 
correlations: 

W.jit -s)^ {V,{t) V,{s)) - {V,{t)) {V^{s)) , z, J - 1, 2, 3 . (2.4) 

Since the generalized hamiltonian V{t) in (2.3) involves stochastic variables, the den- 
sity matrix R(t), solution of the equation of motion (2.1), is also stochastic. Instead, we 



^ Although within a mean-field (MSW) approach, this possibility has recently been 
reconsidered in [15-17] and found compatible with present experimental data. 
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are interested in the behaviour of the reduced density matrix p{t) = {R{t)) which is ob- 
tained by averaging over the noise; it is p{t) that describes the effective evolution of the 
neutrinos in the medium and allows the computation of relevant transition probabilities. 
By making the additional assumption that neutrinos and noise be decoupled aX t — 0, so 
that the initial state is p(0) = (-R(O)) = -R(O), a condition very well satisfied in typical 
situations, an effective master equation for p{t) can be derived by going to the interaction 
representation, where we set: 

R{t) = e**-f^» R{t) e-'*^° , a{t) = e^*^° a 6"^*^° , Lt[ ] = -i\v{t) ■ a{t), 1 . (2.5) 

By averaging R{t) over the noise, we get the reduced density matrix p{t) = {R{t)) in the 
interaction representation: it is convenient to operate on the standard series expansion of 
R{t), so that: 

oo 
k=0 



p(t)=M[p(0)]^ J]ivf)[p(0)] (2.6a) 



(k) 

where the terms A^^^ ' are explicitly given by 

Ni'\pm = / ds^ r ds2--- dsu (L,,L,,---L,J[p(0)] , (2.66) 

Jq Jo Jo 

with N^^^ = 1 the identity, iV^^°^[p(0)] = p(0). The resulting series is a sum over multiple 
integrals of correlators {Lg^Lg^ ■■■Lg^), that is of averages over the noise of successive 
commutators with respect to the stochastic operators V{t) = V{t) ■ a at different times. 
The density matrix p(0) is not averaged over due to the assumption on the initial state 
p(0) = p(0). 

In order to arrive at a more manageable time-evolution, we use a techique [23], of 
which we give a brief account below, that leads to a so-called convolution-less master 
equation and is based on the hypothesis of weak coupling between system and stochastic 
environment. The first step is to write the formal inverse of the map Aft in (2.6a), i.e. 
p(0) =Aff^[p{t)], so that: 

-1 _ /"i ^ atW ^ /\r(2) ^ . . \ _ 1 _ atW _ ^r(2) _ r atW\2 



Aff' =(1+ N^'> + Ni'> + . . .j = 1 - N^^') - Nl^> - {Nl'y , (2.7) 

where only terms containing up to two-point correlation functions have been indicated. 
Further, denoting with N^^^ the time-derivative of N^^^ , it follows that the reduced density 
matrix in the interaction representation satisfies the equation of motion 

^ = ^W'')l = ^^"l«')l = K' + (i*f'-i*.'"iV<") + ...}«*)] . (2.8) 

Since the interaction of the travelling neutrinos with the medium is weak, one can 
focus the attention on the dominant terms of the previous expansion, neglecting all con- 
tributions higher than the second-order ones. Further, since the characteristic decay time 
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of correlations in the medium is by assumption much smaller than the typical time scale 
of the system, the memory effects implicit in (2.8) should not be physically relevant and 
the use of the Markovian approximation justified. This is implemented in practice by 
extending to infinity the upper limit of the integrals appearing in iV(2) and iV(i).[3-5] 
By returning to the Schrodinger representation, one finally obtains [24] 

^ = -i[H,p{t)]+L[p{t)] , (2.9a) 

where 

H = Ho + Hi+H2 = ft-a , (2.96) 

Up] = 2 ^ '^(^ipf^j - , p}] ■ (2.9c) 

The effective hamiltonian in matter, H, differs from the one in vacuum, Hq, by first order 
terms (coming from the piece N^^^ in (2.8)) depending on the noise mean values: 

= {V{t)) ■ a , (2.10) 
and by second-order contributions (coming from the second-order terms in (2.8)), 

3 

H2 = ^ €ijk Cijak , (2-11) 

i,j,k=l 

involving the noise correlations (2.4) through the time-independent combinations: 



Cij 



3 /.Qo 

V/ dt Wik{t) Ukji-t) , (2.12) 



where the 3x3 orthogonal matrix U{t) is defined by the following transformation rule: 
^itHo fj. ^-itHo — Uij{t) aj. On the other hand, the contribution L[p] in (2.9c) is a 

time- independent, trace- preserving linear map involving the symmetric coefficient matrix 
dj = Cij + Cji. It introduces irreversibility, inducing in general dissipation and loss of 
quantum coherence. Altogether, equation (2.9) generates a semigroup of linear maps. 
Ft : p(0) I— > p{t) = /f[p(0)], for which composition is defined only forward in time: 
Ft o Fg = Ft+s, with t,s > 0; this is a very general physical requirement that should be 
satisfied by all Markovian open system dynamics.''' The set of maps Ft is usually referred 
to as a quantum dynamical semigroup. [3-7] 



^ Notice that the procedure of averaging transition probabilities over random matter 
profiles as performed in [25] is not compatible with this basic evolution law, and therefore 
it is a pure phenomenological device, hardly amenable to a rigorous theoretical treatment. 
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The typical observable that is accessible to the experiments is the probability Vv^^v^ (t) 
for having a transition to a neutrino of type at time t, assuming that the neutrino has 
been generated as ai t = 0. In the language of density matrices, it is given by: 



(2.13) 



where Pv^it) is the solution of (2.9) with the initial condition given by the matrix Pv^^) = 
p^^ = \i>f,){i'g\^ while py^ — 1 — py^. By expanding the neutrino density matrix in terms 
of the Pauli matrices and the identity ctq, p = [ctq + p- ^]/2, the linear equation (2.9a) 
reduces to a diflFusion equation for the components pi, p2, ps of the vector p: 



dm 

dt 



= -2np{t) 



(2.14) 



the entries of 3 x 3 matrix H can be expressed in terms of the coefficients fij and Cij 
appearing in the hamiltonian and noise contribution in (2.96), (2.9c): [26] 



n = 



a b + 0,3 c — O2 
c + n2 7 



with a = C22 + C33, a = Cii + C33, 7 = Cn + C22, b = -C12, c = 
solution of (2.14) involves the exponentiation of the matrix H., 



(2.15) 



-Ci3, /3 = -C23. The 



p{t) = M{t) m 



M{t) = e 



-2m 



so that the transition probability in (2.13) can be rewritten as 

3 



1+ E^/^'-. -^^^w 



1 r 
2 



l-M33{t) 



(2.16) 



(2.17) 



Indeed, taking the standard form of the Pauli matrices with respect to the orthonormal 



basis 



and 



I' I,} = 



then p^^ = 



1 + (73 , 1 

and py^ = - 



0-3 



2 ^ 2 

When correlations in the medium are negligible, i.e. the combination in (2.4) are 
vanishingly small, equation (2.9) describes standard (MSW) matter effects, for the presence 
of matter is signaled solely by the shift Hi in the effective hamiltonian. In this case, the 
neutrino-medium interaction is dominated by coherent forward scattering, and, in absence 
of flavor changing effects, the stochastic vector field in (2.3) results oriented along the 
third axis, whence Hi = Aa^, where A = {Vs{t)) = Gpn-e/V^ gives the extra energy 
contribution that electron neutrinos receive when travelling in ordinary matter (Gp is the 
Fermi constant, while He represents the electron number density in the medium). As a 
consequence, the transition probability in (2.17) can be expressed in terms of a modified 
frequency ujm and mixing angle Om in rnatter. 



sin^ 26* + (1 - A/Anf cos^ 29 



1/2 



sin 29 m 



to 



sin 26* , 



(2.18) 
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Afi = a; cos 2^ being the value of A at resonance. In fact, the assumption of neghgible 
correlations amounts to considering in equation (2.14) a matrix Ti. of the form 




n=\-n3 Hi , ni = a;sin(2^) , l]3 = ^-a;cos(26») . (2.19) 



This matrix can be easily exponentiated as in (2.16), 



M{t) 



/ n,+03cos(2u,Mt) _gasin(2a;Mt) ^(1 - cos(2a;Mt)) ^ 
^^(l-cos(2a;Mt)) -gtsin(2c.Mt) "^+"? j(^'^Mt) j 



(2.20) 



whence the explicit form of the element M.33{t) yields the familiar expression: 

P^^^^^ (t) = sin^ 2^M sin^ LOut . (2.21) 

The situation can significantly change for neutrinos immersed in a fluctuatingjnedium; 
while travelling in it, they encoimter matter fluctuations, whose correlations Wij{t — s) 
determine the dissipative contribution in (2.9c). In a typical bath at finite temperature, the 
correlation functions assume an exponentially damped form; therefore, one can generically 
write: 

Wij{t-s) = Wij e-^^^l*-^l , (2.22) 

with Wij and Ay time- independent, real coefiicients, with Ay > 0. Further, as discussed 
before, the stochastic medium fluctuates on time intervals much shorter than the typical 
neutrino "free" evolution time scale 1/u), so that the decay parameters Xij must be much 
larger than the vacuum frequency u). This fact allows neglecting all contributions higher 
than the first-order one in the ratio uj/Xij while evaluating the coefiicients C^- in (2.12). For 
generic correlations as in (2.22), these coefficients, and therefore the entries of the matrix 
Ti, in (2.15), are all nonvanishing. However, the parameters o, b, c, a, /3, 7 describing 
matter decoherence effects are not all free: as we shall see, physical consistency requires 
them to satisfy certain inequalities; in turn, these constraints refiect some fundamental 
characteristics of the matter-neutrino interactions. 

We shall now discuss some interesting cases of the master equation (2.9), corresponding 
to specific physical realizations of the medium through which the neutrinos propagate. 



3. GENERALIZED MSW DYNAMICS 

The simplest instance of a stochastic medium corresponds to ordinary matter with 
density fluctuations, where only the propagation of electron neutrinos is affected. It gen- 
eralizes the familiar MSW mean fleld treatment by adding to it decoherence effects. In 
this case, the stochastic hamiltonian in (2.3) becomes diagonal and, without loss of gen- 
erality, only the stochastic fleld V3{t) can be taken to be non- vanishing; the neutrinos 
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are still forward scattered by the medium, although no longer in a coherent way. This is 
situation discussed in [9-14], where however the density fluctuations in the medium are 
taken to be exactly 5-function correlated. This is a highly idealized assumption, that can 
hardly be reproduced in ordinary conditions. Instead, the much more realistic exponential 
ansatz (2.22) will be used here, where the only nonvanishing correlation strength and decay 
constant are W33 = W and A33 = A, respectively. 

The noise contributions in (2.9) can be explicitly computed; within our approxima- 
tion, one finds that only the entries C23 and C33 of the coefficient matrix in (2.9c) are 
nonvanishing, 

a = C33 = ^, /3 = -C23=^sin2^, (3.1) 

while the hamiltonian contribution Hi is proportional to 0-3 (the standard MSW piece) 
and H2 to CTi: 

l]i =a;| 1 + ^ I sin2^ , 112=0, fla = -co{ 1 ^ \ cos29 . (3.2) 




A2 

Surprisingly, the dynamics generated by (2.9), or equivalently (2.14), with these coefficients 
appears to be physically unacceptable. 

As mentioned at the beginning, any density matrix must be a positive operator [i.e. its 
eigenvalues should be non-negative) in order to represent a physical state: its eigenvalues 
have the physical meaning of probabilities. Therefore, any time evolution needs to preserve 
this property, otherwise an initial state would not be mapped to another state at a later 
time. This is precisely what happens when the neutrino evolution in the medium is modeled 
by (2.9) with dissipative parameters as in (3.1). In fact, the probability V{t) for having a 
transition from an initial neutrino state p(0) to its orthogonal state p± = 1 — p(0) at a later 
time t is given by the first equality in (2.17), with the substitutions pi^^ —>■ p(0), p^^ — > p±. 
Since V is initially zero, its time derivative must be positive at t = 0, otherwise we would 
have physically unacceptable negative transition probabilities as soon as t > 0. A simple 
computation gives: V{0) = Yl'i j=i p{^)i'^ij p{^)j ^ 0, and since this must be true for any 
initial state, physical consistency requires the symmetric part of the matrix Ti in (2.15) to 
be positive. One easily sees that this is impossible with the assignment in (3.1).^^ 

In the case of a 5-correlated medium, the parameter /? identically vanishes and no 
inconsistencies arise; however, as mentioned before, this choice is not supported by strong 
physical motivations and appears just a mathematically convenient simplification. By 
naively relaxing the 5-correlated assumption, one ends up with the simple stochastic system 
discussed above, which turns out to be seriously fiawed. As a consequence, modelling 
matter fluctuations only in terms of electron density is physically untenable and indicates 
that in order to consistently describe neutrino oscillations in random matter more complex 
situations need to be analyzed, involving a richer covariance structure than with V{t) = 
{0,0,V3{t)). 



^ An example of the emergence of negative transition probabilities is explicitly provided 
in the Appendix. 
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Alternatively, instead of the random matter model one may question the approxima- 
tions used in deriving the master equation (2.9), and precisely the weak coupling hypothesis 
and the markovian limit. However, the first assumption appears rather well satisfied in 
the case of the neutrinos, as they interact very weakly with matter, while the markovian 
approximation is justified by the physically motivated choice of rapidly decaying matter 
correlations: A ^ a;. In reality, once a slightly generalized model of random medium is 
adopted, the master equation (2.9) results perfectly adequate to consistently treat deco- 
herence phenomena in neutrino matter oscillations. 



4. DIAGONAL CORRELATIONS 



When the components of the stochastic field V{t) are all nonvanishing, the noise 
hamiltonian in (2.3) is no longer diagonal: in this case, while travelling in the medium, 
all neutrino species undergo incoherent scatterings, in general involving not exclusively 
the forward direction; this may happen only in presence of flavor changing interactions. 
However, as a minimal extension of the previously treated case, we shall assume Vi (t) and 
V2{t) to have zero mean, so that the hamiltonian correction Hi contains only the standard 
MSW contribution, and further take the correlation functions in (2.22) to be diagonal: 



Wij{t-s) = W,e- 



-\i\t- 



(4.1) 



In addition, for simplicity we shall consider situations for which the ratios Wi/\i are all 
equal to a common factor W > 0; in this case, the parameters appearing in (2.15) take the 
form: 



a 



a 



7 = 4W , h = 2ujW 



(---] 



cos 2^ , c = , P = 2ujW 



1+2W 




sin 29 , = , fia = —uj 



A 



A 



+2W 



R 




(4.26) 



and the master equation (2.9) can be exactly integrated. Notice that the request of pos- 
itivity of p{t) for any t > now requires > + P^, condition that is always satisfied 
by the original hypothesis of fast decaying matter correlations: A^ ^ uj. Even more, this 
inequality guarantees not only the positivity of the evolution generated by (2.9), but ac- 
tually a stronger attribute, that of "complete positivity" . [3-5] This property is crucial in 
assuring the consistency of any generalized, dissipative dynamics in all possible physical 
conditions and should always be imposed in place of simple positivity to avoid possible 
inconsistencies in the treatment; [26] it is reassuring that it emerges naturally from our 
simple model of random matter, without the need of further assumptions. 

The transition probability V^^^y^ in (2.17) can be explicitly computed and cast in 
the simple form: 



V.^^.^it) = ^ (l - e"'"*) + e-''^' sin^ 26 sin^ Vtt , 



(4.3) 
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where O = [Of + n| — 6^ — Z?^]-*^/^ is the modified osciUation frequency, while sin^ 29 = 
{Qi — < 1 defines a new mixing angle (notice that the absence of the parameters c 

and O2 is due to the assumptions that led to (4.2a, b)). In comparison with the standard 
result in (2.21), one sees that the presence of a random medium introduces exponential 
damping terms and further modifies the neutrino effective masses and mixing properties; 
a resonance enhancement is still present for A = Ar, but its effectiveness is reduced by 
the damping factors. This is even more dramatic at large times, where the decoherence 
effects dominate: the neutrino state p is driven to the totally mixed state a"o/2 and the 
transition probability approaches its asymptotic 1/2 value. 

These conclusions apply to neutrinos travelling in uniform random media. When 
this is not the case, the neutrino total time evolution results from the composition of 
arbitrarily many partial evolutions corresponding to media with uniform properties, but 
in general of different thicknesses; then, the complete evolution matrix Ai{t) as defined 
in (2.16) will be the result of the composition of the corresponding ones pertaining to the 
various media (a simple example is given in the Appendix). Nevertheless, for slowly varying 
conditions, this composition can be well approximated by its adiabatic expression, obtained 
by the instantaneous diagonalization of the now time-dependent matrix Ti in (2.15) and 
the assumption that the neutrino states evolve as one of its eigenstates.^ Within this 
approximation and neglecting fast oscillating terms, the averaged transition probability 
can be cast in the following form: 



-swt 



1 



n 



1 



A 



4W 



cos 20 



with 



7^ 



4W 
■ Ao 



1 + 



A2 4W 



A3 A 



tan^20+ 1- 



A 



A 



R 



A 4W 



A 



+ 



R 



A, 



1+ 



A2 4W 



Ai Ai 



(4.4) 



1/2 



(4.5) 

With respect to standard, familiar expressions, the action of the stochastic medium 
is signaled by the presence in the second term of a modified weight and a damping factor; 
these additional contributions depend on the ratios of the three matter-correlations decay 
constants A^ and the corresponding strength W. Although in the weak-coupling regime 
one expects W <ti u, the decay constant ratios need not be small. Therefore, the behaviour 
of (4.4) as a function of the neutrino energy can sensibly differ from the one obtained in 
absence of decoherence effects (concrete examples are shown in Fig.l). 

Of particular interest is the application of (4.4) to the solar neutrino case, where 
Am^ and 9 can be taken to assume the best fit values obtained in recent data analysis 
{e.g. see [27, 28] and references therein); thanks to the availability of a larger decoherence 
parameter space, the electron surviving probability Vv^-^v^ = 1 — Vv^^vy, is found to differ 
not only from the standard, noiseless expression, but also from those obtained with 8- 
correlated fluctuating matter as reported in [13, 14] (for a comparison, see Fig.l). These 



't' Possible hoppings among the instantaneous eigenstates can also be easily included; 
for simplicity, we ignore them here. 
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results, together with the still present uncertainties in the fluctuating behaviour of the 
solar matter, appear to open concrete possibilities for an experimental study of matter 
induced effects in neutrino oscillations that go beyond the standard MSW phenomenology. 



5. DISCUSSION 

In the most general situation, the correlations in the stochastic medium have the 
form (2.22) and thus all the entries of the matrix in (2.12) result nonvanishing; as a 
consequence, all second order pieces in the effective hamiltonian (2.96) as well as in the 
dissipative part (2.9c) will contribute to the master equation (2.9o). Further, the first 
order mean field approximation in (2.10) will no longer be diagonal, taking into account 
the presence of possible fiavor changing interactions. [15-17] 

Nevertheless, even in this very general case, the corresponding matrix Ti. in (2.15) can 
not result totally generic: as already pointed out, the positivity of the evolved state p{t) 
must be preserved under all circumstances; this is guaranteed by the mentioned condition 
of complete positivity of the evolution generated by (2.9). This property requires the 
positivity of the matrix in (2.9c) and as a consequence imposes certain inequalities 
among the dissipative parameters in (2.15) (see [29, 30] for explicit expressions). These 
conditions are certainly of help in restricting the parameter space needed to describe a 
totally generic random medium. 

Even with these constraints, no simple, exact analytic expressions for the transition 
probability 'Pu^^v^{t) in (2.17) can in general be given. However, as discussed before, 
second-order matter contributions to Ti. are small with respect to the vacuum frequency 
a;; therefore, in solving (2.14) one can integrate the hamiltonian dynamics exactly, while 
treating the dependence on a, 6, c, a, /3, 7 in perturbation theory. ^ In this way manageable, 
approximate expressions for the transition probabilities can be obtained. Having now 
at disposal a larger parameter space, their form involves multiple damping factors and 
oscillation phases, showing possible larger deviations from the standard behaviour. In 
addition, notice that in order to describe neutrino mixing in a generic random medium 
two mixing angles, 6 and are in general needed: they parametrize the components of 
the unit vector 0^/0 = ( cos sin 2^, sin sin 2^, — cos2^), with VL — which identifies 

the effective hamiltonian H — Vt ■ a 'm (2.96). 

Actually, in presence of Majorana neutrinos, also in vacuum the most general mixing 
matrix involves two angles, 9 and so that the explicit expression of the free effective 
hamiltonian Hq in terms of these angles is as for H above. Although for oscillations in 
vacuum involving only two species of neutrinos the angle disappears from all observables, 
this is no longer true in presence of matter induced decoherence effects. Indeed, one can 
directly check that the transition probabilities explicitly depend on unless the dissipative 
parameters a, 6, c, a, 7 are all zero; at least in principle, it is therefore possible to 
distinguish between Dirac and Majorana neutrinos by studying their oscillations in random 

^ The effects of the hamiltonian corrections to the free motion are in general not small, 
in particular near resonance; this is why no approximation is allowed in the evolution 
generated by the effective hamiltonian (2.96). 
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matter. The detailed analysis of such of dependence is certainly beyond the scope of the 
present investigation and thus, in order to keep the treatment as simple as possible, in the 
previous discussions wc have tacitly assumed the neutrinos travelling in matter to be of 
Dirac type, setting (p — from the beginning. 

As a final remark, let us mention that master equations of the type (2.9) generate the 
most general open system dynamics compatible with a semigroup composition law and the 
requirement of complete positivity, and as such can be applied to model in a physically 
consistent way a wide range of phenomena. [3-7] In particular, they have been recently used 
in order to describe dissipative effects induced at low energies by the dynamics of funda- 
mental objects (strings and branes) at a very high scale, typically the Planck mass. [29, 30] 
These string induced decoherence effects may modify the pattern of neutrino oscillations, 
and in principle interfere with the phenomena described above. Nevertheless, besides be- 
ing very small, they affect in equal manner all types of neutrinos, so that they can be 
isolated from the matter-induced effects by analyzing data taken in different experimental 
conditions. 



APPENDIX 

In order to show that negative probabilities arise in experimental accessible obscrvables 
once the naive model of density fluctuating matter discussed in the text is adopted, one 
needs to combine neutrino propagation in vacuum with that in the medium. Consider a 
neutrino, created as Ue, that propagates for time t' in vacuum, then enters the random 
medium in which stays for a time t, and is finally detected after having travelled again in 
vacuum for a further time t". The probability Vy^^u (r) of finding a neutrino of type 
at the final time r = t' + t + t" can be expressed as in (2.17), where the total transition 
matrix M{t) is now the product of three terms, M{r) = Mo{t') ■ M{t) ■ Mo{t"), the 
middle representing the propagation in the medium with parameters as in (3.1) and (3.2), 
while the outer two the "free" motion in vacuum, generated by the hamiltonian (2.2). 

When the vacuum evolution time t' is chosen to be very short, such that sin ujt' = 
13/ (2a sin 2^) = uj/A\, the state of the neutrino entering the medium is p_ = [ctq — p_ ■'?]/2, 
where p- = —Aio{t')p,^^{0) coincides with the eigenvector of the dissipative part of Ti. 
relative to its negative eigenvalue. Similarly, with the same choice also for t" , one finds 
•^o{—t")Piyi_,{'^) = P-) so that when exiting the medium the neutrino is found in the state 

p+ = 1 — p_, orthogonal to p_. With these conditions, one has: Vi,^-h,^{t) = V >.+ (t), 

and near resonance, one explicitly finds: 



-at 



^ sin20t 
cos 2nt + D — — 



where D — [a^/4-l-/3^]^/^ and Q, = [a;^ — D^]-*^/^; this expression indeed assumes unphysical 
negative values for sufficiently small times: V,^^->,^^{t) ~ (a/2 — D)t. 
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0.1 1 ^° A/Ap, 



Figure 1. Behaviour of electron neutrino mean survival probability V,y^^,^^ as 
a function of the neutrino energy (through the ratio A/Aji), for sin^ 2^ ~ 0.8, 
density dominated matter fluctuations, Ai, A2 S> A3, and different correlation 
strengths, W/A2 — 10"'^ 10""^. The lower starting (black) curve corresponds to 
the case of noisless matter (standard MSW effect), while the remaining (colored) 
ones show the effect of the stochastic fluctuations. The initial gap among the 
group of curves is due to the presence of the decoherence driven damping factor. 
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